OPTIMAL ERROR ESTIMATES AND ENERGY CONSERVATION 
IDENTITIES OF THE ADI-FDTD SCHEME ON STAGGERED GRIDS 
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Abstract. This paper is concerned with the optimal error estimates and energy conservation 
properties of the alternating direction implicit finite-difference time-domain (ADI-FDTD) method 
which is a popular scheme for solving the 3D Maxwell equations. Precisely, for the case with a 
perfectly electric conducting (PEC) boundary condition we establish the optimal second-order error 
estimates in both space and time in the discrete _ff ^-norm for the ADI-FDTD scheme and prove the 
approximate divergence preserving property that if the divergence of the initial electric and magnetic 
fields are zero then the discrete L'^-norm of the discrete divergence of the ADI-FDTD solution is 
approximately zero with the second-order accuracy in both space and time. A key ingredient is 
two new discrete energy norms which are second-order in time perturbations of two new energy 
conservation laws for the Maxwell equations introduced in this paper. Furthermore, we prove that, 
in addition to two known discrete energy identities which are second-order in time perturbations of 
two known energy conservation laws, the ADI-FDTD scheme also satisfies two new discrete energy 
identities which are second-order in time perturbations of the two new energy conservation laws. 
This means that the ADI-FDTD scheme is unconditionally stable under the four discrete energy 
norms. Experimental results are presented which confirm the theoretical results. 
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1. Introduction. The alternating direction implicit finite-difference time- 
domain (ADI-FDTD) method was first proposed in 2000 in [121 [55] for the three- 
dimensional Maxwell equations with an isotropic medium. This method is based on 
the staggered grids of Yee's scheme [27] and splitting of Maxwell's equations, consists 
of only two stages for each time step and is unconditionally stable. The ADI-FDTD 
scheme is a popular scheme for solving the three-dimensional Maxwell equations and 
is applicable to a wide range of problems in computational electromagnetics (see [25]). 
So far, the ADI-FDTD scheme has been intensively studied by many authors (see, 
e.g. [21 [3 H H H [3 [m [201 m [24, 26, 28, 30i), such as the stability and dispersion 
analysis using Fourier analysis ^ |8, |20, i23, |251 [30] , the unconditional stability and 
convergence analysis by the energy method [5] [5], its accuracy and efficient algorithms 
[31 [3 mi US] and its applicability with perfectly matched layer absorbing boundary 
conditions [13] ■ Meanwhile, motivated by the success of the ADI-FDTD scheme, 
there is also considerable interest in developing new unconditionally stable schemes 
for solving the Maxwell equations, based on splitting of the Maxwell equations (see, 
e.g. [Il[2l[5l[l0l[ni[l2l[l5l[l6l[22]), where the stability and convergence of these 
schemes were studied in the discrete LP' norm. 

However, optimal error estimates and energy conservation and divergence preserv- 
ing properties of the ADI-FDTD scheme have not yet been studied although these 
properties are vital for the successful application in practical problems of the ADI- 
FDTD scheme. The purpose of the present paper is to conduct such a study. In 
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particular, for the case of perfectly conducting (PEC) boundary conditions we will 
establish optimal error estimates in the discrete H^-noTm for the ADI-FDTD scheme 
and prove that the ADI-FDTD scheme satisfies four discrete energy identities which 
are second-order in time perturbations of four energy conservation laws, employing 
the energy method. We also prove an approximate divergence preserving property 
with a second-order accuracy in both space and time for the ADI-FDTD scheme. In 
doing so, a crucial tool is two new discrete energy norms (cf. Theorems 13 . 21 and 13 .31 in 
Subsection I3.3P which are second-order in time perturbations of two new energy con- 
servation laws for the Maxwell equations (see Theorems 12 . 1 1 and 12 . 2l below') introduced 
in the present paper. 

The remaining part of the paper is organized as follows. In Section [21 we derive 
two new energy conservation laws for the Maxwell equations with the PEC bound- 
ary conditions and introduce two well-known energy conservation laws. In Section 
[3] we recall the ADI-FDTD scheme, define some new energy norms and prove that 
the ADI-FDTD scheme satisfies four energy identities which are second-order in time 
perturbations of the four energy conservation laws for the Maxwell equations intro- 
duced in Section [21 This means that the ADI-FDTD scheme is unconditionally stable 
under the four discrete energy norms. By making use of the two new discrete energy 
norms we derive optimal second-order error estimates in both space and time under 
the discrete norm for the ADI-FDTD scheme (see Theorems 14. 2[ 14.31 and 14. 5p in 
Section [4l and the approximate divergence preserving property (see Theorem 15. 1|) in 
Section [5l Note that it is a superconvergence result that the ADI-FDTD scheme is 
second-order convergent in space under the discrete norm. Experimental results 
are presented in Section [HI which confirm the theoretical results. 

2. Maxwell's equations and their conservation properties. In this section 
we introduce four conservation properties of Maxwell's equations, of which two are 
new. These conservation properties will be used to study the ADI-FDTD scheme in 
this paper. 

2.1. Maxwell's equations. Maxwell's equations governing the propagation of 
electromagnetic waves can be written as a system of partial differential equations (see 

my- 

(2.1) VxE = - — , VxH = J + — , V-B = 0, V-D = p, 

where E is the electric field, H is the magnetic field, B — /iH is the magnetic flux 
density, D = eE is the electric displacement, J = crE is the current density, e is the 
electric permittivity, /i is the magnetic permeability, a is the electric conductivity and 
p is the volume density of electric charge. 

In this paper we consider Maxwell's equations in a lossless and isotropic medium 
without a source field. In this case, /x and e are constant and p = ct = 0, so Maxwell's 
equations become 



(2.2) 
(2.3) 
(2.4) 
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where E = {Er^, Ey,Ez) and H = {Hr^, Hy, H^) denote the the electric and magnetic 
fields, respectively. For simplicity we consider the perfectly electric conducting (PEC) 
condition on the boundary dil of the cubic domain il = [0, 1] x [0, 1] x [0, 1]: 



(2.5) 



n X E = 0, n • H = 0, on (0, T] x dn, 



where T > and ft is the outward normal vector of dil. We also need the initial 
conditions: 



(2.6) 



E(0, a;, y,z) Eo(a;,2/, z) and H(0, x, y, z) = Ho(a;, y, z) 



It is well known that, for suitably smooth data, the problem 
solution for all time (see [17]). 



has a unique 



2.2. Energy conservation laws in a lossless medium. When the medium 
is lossless and isotropic, we have the following two new laws of energy conservation. 

Theorem 2.1. (Energy Conservation I) Let E and H he the solution of the 
Maxwell equations \2.2^ — in a lossless and isotropic medium satisfying the PEC 

boundary conditions \2.5\ . Then for w — x,y, z, 



(2.7) 





9E 


2 


OH 






dw 


+ /i 


dw 





dxdydz = Constant. 



Proof. We only prove (|2.7p for the case with w = x. The cases with w = y, z can 
be shown similarly. 

Differentiating the first two equations in (j2.ip with respect to x, we have 



V X 



SE 

dx 



dxdt ' dx dxdt 



Taking the scalar product of the above two equations with dH/dx and dE/dx, re- 
spectively, integrating the resulting equations over fl and adding them together, we 
obtain on integrating by parts that 



(2.9) 



d^inB.) an 92 (eE) dE 
dxdt dx dxdt dx 



f dE da 

dxdydz — ~ {n x — — ) • ——ds. 

Jdn 



Now, from the PEC boundary condition (|2.5p we have n x dE/dx = for y = 0, 1 
or for z = 0, 1. 

Then SIM becomes 



(2.10) 
where 



ld_ 

2dt 



dU 



dx 



dE 



dx 



dxdydz 



To+Ti, 



T, 



dEydH^ dE^dHy 



dx dx 



dx dx 



dydz, i 0, 1. 



Using the Maxwell's equations and the PEC boundary condition it can be shown that 
Tq = Ti = 0. Using the second equation in (|2.2p and the first equation in (j2.3|) we 
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have 

To - 



Jo 



dH^ dEy dH^ dE, 



dz dx 

dH, dEy 
dz dx dy dx 



dy dx 
dH^ dE, 



/dEy dEy dE.dE, 



\ dt dx 
dydz 



dt dx 



dydz 



x=0 



since, by the PEC boundary condition n x E = on dil in ()2.5p . we have dEy/dt = 
dEz/dt = for a: = and for alH > 0, y G fJ. By integration by parts it then foUows 
that 
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0, 



where use has been made of the fact that = for a; = and for all t > 0, y G 
in view of the PEC boundary condition n • H = on dVl in (|2.5p . 

Similarly, it can be shown that Ti — 0. Thus ()2.10p implies (|2.7p with w = x. 
The theorem is thus proved. □ 

It is easy to see that the above proof of Theorem 12.11 still works if E and H are 
replaced by dE/dt and dH/dt, respectively, so we have the following result. 

Theorem 2.2. (Energy conservation II) Let E and H be the solution of the 
Maxwell's equations i2.^ )- [2^ in a lossless and isotropic medium satisfying the PEC 
boundary condition 1^2.5} . Then for w — x,y, z, 



(2.11) 





d^E 


2 




d^H 






dwdt 




dwdt 





dxdydz = Constant. 



As far as we know, the two energy conservation laws (12.71) and (|2.1ip above are 
new. However, the following two energy conservation laws are already known and 
were proved in PO; in particular, the first one is the well-known Poynting theorem and 
can be found in many classical physical books. 

Lemma 2.3. (Energy conservation III) Let E and H be the solution of the 
Maxwell's equations i2.2\) - (2^ in a lossless medium satisfying the PEC boundary 
condition i2.5\). Then 



(2.12) 



/ (ejEp + At|Hp) dxdydz = Constant 
Jn 



Lemma 2.4. (Energy conservation IV) Let E and H be the solution of the 
Maxwell's equations ll2.2\) - (2^^ in a lossless medium satisfying the PEC boundary 
condition i2.5\). Then 



(2.13) 




(9E 



dt 



dU 



dt 



dxdydz — Constant 



We now have four conservation laws (|2.7p , ()2.1ip , ()2.12|) and (|2.13|) for the prop- 
agation of electromagnetic waves in a lossless and isotropic medium. We will prove 
in the next section that the ADI-FDTD scheme keeps these properties with a second- 
order in time perturbation term. 
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3. Energy identities of the ADI-FDTD scheme in H^. Before establishing 
the energy identities for the ADI-FDTD scheme, we first recall the ADI-FDTD scheme 
(see |28j ) and define some discrete energy norms. 

3.1. The ADI-FDTD scheme. The ADI-FDTD scheme makes use of Yee's 
spatial staggered grids (see [27j, which are denoted by fl'^ and defined as: 

n'^ -. {{xa,yi3,z^) I Xa = aAx, yp = /3Ay, = •jAz, a = i,i + 1/2, 

i = 0, 1 • • • , / - 1; /? - J, J + 1/2, J = 0, • • • , J - 1; 7 - fc, fc + 1/2, 
k ^ 0, ■ ■ ■ , K - 1; xo = yo = zq ^ 0, XI = yj = ZK ^ 1}, 

where Ax, Ay and Az are the mesh sizes along the x, y and z directions, respectively. 
In the ADI-FDTD scheme, the field components are defined on different subsets of 
fl'^ (see the subscripts in the equations of this scheme below). 

For a positive integer N let At = T/N be the time step and define =: nAt 
with n = 0,1,---,A". Denote by u™^ ^ the grid function defined on the point 
(rn At, aAx, /3 Ay, 7A2;). Then define 



7 ^a,/3,7 



)/At, 



^v'^a.,i3,'y ^^'^ '^^"™/3,7 defined similarly. Denote by E^^j; ,^ and H^a.,/^,^ 

the approximations of the electric field E^it"^ , Xa,yp, Zj) and the magnetic field 
Hu,{t"^,Xa,yi3,Zj), w = x,y,z. Then the ADI-FDTD scheme (see [28]) is defined 
in the following two stages. 
Stage 1: 



(3.1) 
(3.2) 
(3.3) 



rpn Tpn 

Ex - Ex 




At 




E^ - E: 




At 
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SyH" — 5zHy 
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rrn rrn 
Ex - tlx 
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At 
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m - h: 




At 


^ 2^ 



SzHx — SxH" 



i,3,k 



whereafter I :— I + 1/2 for I — n, k and F\a^i3^-y means each term of the expression 
or the equation F has the subscripts a, /3, 7. 
Stage 2: 
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The PEC boundary condition (12. 5p can be discretized as: 
(3.7) 
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where m = n or n + 1/2 denotes the time levels and k are integers in their valid 
ranges. We also need the discrete initial conditions which are obtained by imposing 
the initial condition (12.61) at i = 0: 



EO^^_^ = Eo(aAa;,/3Ay,7Az), 



Ho 



Ho(aAa:,/3Ay,7Az). 



3.2. Discrete energy norms and notations. We will need some dis- 



crete energy norms and notations. For U = {Uxj, 



(^a;ij,fc,K;j V^ij^fc), define the discrete norms || • and |j • \\h as: 



null 



I-l J-1 K-1 



I-l J-1 if-1 



J-1 J-1 K-1 



p(^)uk,^. + E E E ^(u)f/^,,,. + E E E ^wt/^..,. 

i = j = l fc = l 1=1 j = fc = l i = l j = l fc=0 

E E E p(y)v^.u + E E E ^(v)^,,.. + E E E pi^whu 



2=0 j=0 fc=0 



i=0 j=0 fc=0 



i=0 j=0 fc=0 



Aw 



Av 



where Av = AxAyAz, P(E) — e, the electric permittivity, and P(H.) ~ the 
magnetic permeability. For the electric field E and the magnetic field H we also need 
the following norms near the boundary grids: 



J-lA'-l J-1A--1 

EE E ^^^L,,, + EE E 



j=0 fc=l i' = ij-i 
J-1 K-1 



AyAz, 



j=0 fe=0 



\i,j.k 



\l-l,3.k 



AyAz. 



The norms 1 1 E 1 1 j , 1 1 E 1 1 ^ , 1 1 H 1 1 j and 1 1 H 1 1 are similarly defined by changing the indices 
j, fc, the factor 1/Ax and the area element AyAz into the corresponding ones. 

Furthermore, for simplicity in notations, we introduce some difference operators. 
For a vector-valued grid function U = (Uxa,i3,'y, Uy^ ^ ^, Uza,0,^) with a = i or i + 1/2, 
= j or j -I- 1/2, 7 = A: or fc -I- 1/2 define 

(5.u,U = {&wUx,5wUy,5wUz)\a,l3,'y, w = x,y,z, 

SlV^{6yUz,SzUx,SxUy)\ Sl^U=(SzUy,SxUz,SyUx)\ a,/3,7 ■ 

By composition of operators it can be easily derived that 

(5l (5''U = {Sy6xUy,Sz6yUz,SxSzUx)\a,l3,-f, 
S262IJ = i5zSxUz,SxSyUx, SySzUy)\a,li,J. 

3.3. Four discrete energy identities for the ADI-FDTD scheme. We now 

derive four discrete energy identities for the ADI-FDTD scheme, which are second- 
order in time perturbations of the four energy conservation laws for the Maxwell 
equations. To this end, we need the following result on summation by parts. 

Lemma 3.1. (Summation by Parts) Let M > 1 be an integer. For any sequence 
{Wra+i/2}mZl let S is a difference operator defined by 5Wm = {W,n+i/2-Wm-i/2)/h, 
where h is a positive number. Then, for any two sequences {t^m+i/2}m=o '^"■'^ 
{ym}m=o ™e have 



M-l 



M-1 



(3.; 



E f^m+l/2'5^m+l/2 — C^M-1/2^M " t^l/2^0 " E ^^SUm- 



m=0 



m— 1 
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Theorem 3.2. (Energy identity I) Let n > and let E" = 
(EhJ,k,El\J,k'E^^,J,k)' H" = be the solution of the ADI- 

FDTD scheme — \3.6^ with the boundary condition \3.'1\ . Then we have the 

following energy identities: 



|(5^E"+i||| + ||,5„H 



n+l,|2 , i^tf 

\h 



+I|e"-*-^IIlw + I|h"+^|'^ 



(ll^2H 



/iTTn+l||2 



n+l||2 

\l(w) 



h-vTtii+l |i2 



|(5fE' 



rn\\2 



= ||5^E"|li + ||<5^H"||^ 



(3.9) +||E"|'2 



L{w) 



|H" 



\L(w) 



5^(5^E«||2, + p„52^H"|l| 



I ^/ijjn ||2 I ||^/iT7'^'i||2 



\L(w) 



where w — x,y, z, L(a;) = /, L{y) = J, L{z) — K . 

Proof. We only prove p.9p for the case with w = x. The other cases with w = y,z 
can be proved similarly. 

Applying the difference operator Sx to each equation in Stages 1 of the ADI- 
FDTD scheme, rearranging the terms according to the time levels and squaring both 
sides of the six equations thus obtained, we have 



(3.10) e{S,E:+^f + i^(SxSyH:^^f - AtSxE^+hjyH:^'' 



n\2 



{Atl 

As 



(3.11) sis.E;^^y + i^{Sx6,H:+^y ~ AtSxE;;+hj,H:^^' 



{Aty 



n\2 



(3.12) e{5.,E'r~^) 



4e 



[SJ^H^Y ~ AtS^E'^SJ^Hl 
(SJ^Hy^^Y - AtS,E"+hj,Hy' 



e{5xE\ 



n\2 



(Atf 



{SJyH'iY ~ AtS^E^dJyHl' 



(3.13) fiiS^H" 



+ ^-^{5J.E';+'-f - AtSxH:+hj,Ey 



(3.14) fi{SxH;^^) 



lx{5xH:f + IM!(45^i5n)2 _ AtS^H^SJyE^ 

' ^^*^\s.JxErir - AtSxH^^hj.xE:+^ 



2 'J^ 2 ^"-^ 2 



4m 



Afi 



{5J.E': f - At5xH^5J,El 



[Aty 



(3.15) ti{5xHT-f + i-A(^,5^£;r5)2 _ AtS^Hr^SJyE.. 

4^ 



i,3,k+h 



n\2 



4/. 



{S.6,KY - AtS^H^dJ^El 



8 



L. GAO, B. ZHANG 



From the definition of 6x and tlie conditions (|3.7p it follows that 

where m — n, n + 1/2, j' — 0, J, fc' = 0, if. By these boundary conditions and 
summation by parts (Lemma 13. ip . we can see that the sum of all the mixed-product 
terms such as dxE^^S^SuHy {u,v — x,y,z) on the left- and right-hand sides of p.lOl) - 
p.lSp will vanish, respectively, due to cancelation. For example, consider the sum of 
the mixed-product terms on the left hand sides of (|3.15l) . (|3.13p and (|3.14l) over i, j 
and k in their valid range. We have 



/-I K-l 



1=1 j=0 fc=l 
1-2 J-1 X-1 

SxE^SxS^Hy 

i = l i = l 1 



2 ••'^2 •"■^2 



I-l .7-1 K-1 

i=l j=0 1 i,j+5,fc i=l j = l /.— I 

I-2J-1K~1 /-2.7-lX-l 

EEE^^-f^-'^-'^-^^^ 

i=l j=0 fe=0 
I-l J-1 K-1 

SxHyS^s^E" 

1=1 j = l fe=0 i,j,k + ^ 

1 r n 



j=l fe=0 



The sum of the above three terms can cancel the sums of the three mixed-product 
terms on the left-hand side of p.lOp . p. lip and (I3.12p . respectively. Thus, sum up 
each of the six equations above over the valid ranges of their subscripts j,k and add 
the updated six equations together to deduce that 



n l|2 
E 



4/i£ 



.(||5.5rE"||lf + ||5.52H"||^) 



(3.17) 



Aa; ^ ^ 

j = l fc=0 

At 



/-lj,fc+i 



+ Ez 2 d^-ffy 



j = *; = ! 



/-l,3+i,fc 



i,i+i fc- 



AyAz 



where we have used the boundary conditions (13. 7p applied by the operator dx to get 
the last two terms on the right-hand side of the above equation. 

Similarly, from the equations in Stage 2 applied by the operator 5^ we can derive 

that 



n + l||2 , (At)2 



2 



At 
Aa; 

At 



4/ie 



j = l fe=0 

J-1 K-l 



/-lj,fc+i 



ij,fc+i 



(3-18) +£EE(^^""^^^^^" 



j=0 fe=l 



+ Ey"+^5xH,"+^ 



I-lJ+i fc 



j AyAa 
AyAz. 
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We now consider the mixed-product terms on the right-hand side of p.l7p and 
(|3.18p . Taking i — 1 and / — 1 in the first equation in (13.21) and (13.51) we have 

n+i _ At „+i _ p n _ rj n 



E. 



At 



At, 



where i' — 1, / — f . Multiplying both sides of the above two equations by ^/e, squaring 
them and adding the resulting equations give 



2AtE"^^5^Hy^^ = eiE'i+^f - s{E^f 
(3.19) 



(At) 



2 r 



+At( E^+^6yH'^'+^ + E'^SyH, 



i' = 1, 1 -1. 



Similarly, it follows from the second equation in p.ip and p.4p with z = 1, / — 1 that 

At(^E;SxH': + s;+i<5.i/r+i) = £(S;)' - £(s;'+')' + 2AtEy^hM'^^ 
(3.20) _ f^iS^H^y) 

Combining (jXTTI . (|3J8)) . ((3l9l) and (13:20)) gives 



1, /- 1. 



"+i||| + ||<5,H 



■n.+ l \ \2 I rh-|7Ui+l||2 I ||r r/iTT)i+l|i2 



|Mi'E"+^||^ + |15,<5^^H"+^||i) 



i' = l,/-l j = l fc=0 
J-1 A'-l 



AyAz 



+i.E EE[-(^/"^)^ + ^M(w"^)^ 



i' = l,/-l j=0 fc=l 



i' J + 



AyAz 



n||2 I iij jjri||2 , /||^ x''T7i»i||2 i lljf x'>Tjn||2 



ff-^VIlMfE"||^^ + ||M2"H"||i) 



J-l K-l 



-i^ E EE^(^^")^ + ^M(W^^ 



i' = l,/-l j=l k=Q 



(3-21) +^ E EE[+-(V')^ 

i' = l,/-l j=0 fc=l 

where the term i? is given by 

.7-1 if-l 



4/ie 

4^£ 



AyAz 



AyAz + i?, 



E EEK'"^W^ + ^^"W 



^ E E E ^/"^^^^^""^ 



' = !,/-! j = l fc=0 

J-l X-l 



Ax 



i' = l,/-l i=0 fe=l 



Aj/Az. 



AyAz 
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Arguing similarly as in deriving p.l9p we deduce from the second equation in 
and (|3.5p with i ~ i' that 



i' = l,/-l j=0 k=0 ^ 



2 



(3.22) -{t^{H."+'r + ^-j^ei6yE^-+'r} 



AyAz. 



Substituting p.22p into p.2ip gives the identity p.9p with w ~ x, noting the definition 
of the norms ||E||7 and ||H||/. The proof is thus complete. □ 

Note that the proof of Theorem 13 . 2 1 does not depend on the time levels. Thus, if 
we apply the operators St and S^St with w — x,y, z to the equations in the ADI-FDTD 
scheme and repeat the above argument, then we can obtain the following result. 

Theorem 3.3. (Energy identities II) Let n>\ and let E", H" he the solution 
of the ADI-FDTD scheme. Then 



4^£ 

+ ||5,E"+^||i(„) + II W+^lli^ + ^(||<52^<5*H"+^||i(„) + ||<5i%E«+^||i(„; 



= ||(5„,5tE"-^||| + ||5^,5tH"-5||2^ + M!(p„j^j,E"-^||2^ + 
+ ||^,E"-^||i(^) + ||<5,H"-^||i(^) + M!(||4M,H"-i||i(„,) + ||5i'^,E«-^||i(^; 

where w = x,y, z, L(x) = I , L{y) = J and L{z) — K . 

Theorems 13. 21 and 13. 31 are reduced to the following two results, respectively, when 
5x — I (identity operator). 

Theorem 3.4. (Energy identity III) Let n > and let E", H" he the solution of 
the ADI-FDTD scheme. Then 



|E"+i||?. + ||H"+i||2 



(At)! 

E"||| + IIH-lll, + M!(||5^^H"||| + ||5^E"|||, 



^ ^ ^ __(||^2^H"+i||| + ||<S?E"+i|||, 



Theorem 3.5. (Energy identity IV) Let n > 1 and let E", H" he the solution of 
the ADI-FDTD scheme. Then the ADI-FDTD scheme satisfies the following identity: 

||<5*E"+^||| + ||5,H"+^||1, + M!(||5^^5,H"+i||| + II^^J'^^^^^ 

= ||5tE"-3||| + II W-^lll^ + ^(||<5^MtH"-^||| + ||^JW-5||2^). 

Corollary 3.6. The ADLFDTD scheme ([iJ]) - (HJ) with the PEC boundary 
condition \3. 7p is unconditionally stable under the new defined discrete energy norms 
and under the discrete norm. 
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4. Optimal error estimates for the ADI-FDTD scheme. In this section, 
we derive optimal error estimates for the ADI-FDTD scheme under the discrete energy 
norms. We need the following well-known result called the discrete Gronwall's lemma 
(see [14]). 

Lemma 4.1. //{a„}, cind {c„} are three positive sequences with {c„} being 
monotone such that a„ + bn < c„ -|- A X]"=o with A > and oq + &o ^ Cq then 
o-n + bn < c„ exp(nA) for all n > 0. 

For n > 0, w — x,y,z and for a — i,i, l3 = j — k,k let /i -, ~ 

W — where e = (e^^jByje^), h = {hx,hy,hz) is the exact solution to 

the Maxwell equations (|2.2I) - (|2.4I) with the boundary condition (|2.5p and the initial 
condition m and E" = i??,,,, fc), H" = i?;^,^,^, is 

the solution of the ADI-FDTD scheme. 

Theorem 4.2. Let the exact solution e, h satisfy that e e C{{0,T];C'^{Q)) f] 

c\{o,T];C\n)) n c'iio,T];C\n)) n cmo,T];C{n)), h e C{{0,T];C\n)) n 

C\{0,T];C^{n)) n C'^{{0,T];C\n)) n C^iiO,T];C{n)). Then for any n > and 
any fixed T > there is constant C independent of At, Ax, Ay, Az such that 

(4.1) Wd^s'^+'Wl + \Kn"+'fH + ^(II^»'5^^"+'IIh + Wsy^n^'+'Wl) 

4/ie 



< C{{Atr + (Ax)4 + {Ayr + (Az)n, 



where w = x,y, z, L(x) = /, L{y) = J and L{z) = K . 

Proof. We only prove (|4.1I) for the case with w = x. The other cases with w = y, z 
can be proved similarly. 

In order to derive the error equations for the ADI-FDTD scheme and the (5^;— ADI- 
FDTD scheme (i.e. the ADI-FDTD scheme applied by Sx), we need two discrete 
forms of the Maxwell equations corresponding to the two schemes. Denote by e^*, 

{w = x,y,z) the intermediate variables defined by 



"xe^ — 


\5x{e-+' 


+ eL') + 






SxC"^* = 


^^(er^ 


+ ^y) + 


^SxSxih:+' 




SxeT = 


\5x{eT' 


+ e^) + 


^SxSy{h:+^ 




5xK* = 


\5x{h^^' 




—S S ('e"+i 
^^OxOy^e^ 








+ K) 4 


^....(er^ 


^x ) li ,j,fc+i I 


SxK* = 




+ K) + 


—5xSx{e'i+' 





Then by a direct calculation we derive that 

Stage 1 of the discrete form of the Maxwell equations: 

(4.2) ^ = - [^y^^K* - SJxh;) + 5,131 
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5x6^ - 5,e^ _ 1 



M/2 

- Sxe'l 1 



SxSxK)+Sx(3;^~' 

6xhT - SxK ^ 1 _ ^ ^^^.+1 



At/2 jjL 

6xhT - ^ 1 
At/2 /i 



(j^5,er - 5x6x6^) + 



»+5.i+5.fc+5' 



Stage 2 of the discrete form of the Maxwell equations: 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



At/2 

"xCy Ox 

At/2 e 



SJxh?+' 



At/2 
At72 

(4.12) -"^ 



AV2 



(4.13) 



5x/i^+i - SxhT 



At/ 2 



= - (^5ySx 

^{d.SxK 

[SxSxh;* - 6y5xK+^ 

^ ('^z'^xe^* - 5y5xe 

- i^xS. 



"* _ A A /,"+! 



n+1 



A A 



-'X'^X "Xl^X'^y 



+ 6xP2^^ 



+ 5xP: 



n+i 



n+i 



+ 5xi 



+ 5xC'' 



+ 5xi. 



Here, in Stage 1 and Stage 2 of the discrete form of the Maxweh equations Sxfiu 

and 5x£^ ^ {w = x, y, z) arc the truncation error terms of the equivalent form of the 
(5x-ADI-FDTD scheme (sec below). For example, 

{Atf 



Sxi3:^~' 




1 

— — 






£ 




= SxSthx^^ 


1 









A "+2 S 



{Atf 

AjJLS 



SxSySxStey^ '' 



5x5z5x5th/l^' 



The other terms can be obtained similarly. Then the equivalent form of the Ja,-ADI- 
FDTD scheme can be easily derived by combining all the equations in the ^j,-ADI- 
FDTD scheme and canceling the field values at the intermediate time levels t"+2 as 
follows: 



5x5tEx^^ 


= -Sx^t 

e 






+ ^^^Sx5y5x6tE;+^ 


i,3,k 


5x5tEy^^ 




SzHx^'^ 


- SxHz^-' 


+ ^^SxSzSyStE^^-- 




SxStE^^'' 


= -SxSt 

e 


A 17"+ 5 
UxJ^y 


- SyH:+--' 


+ ^-^5x5x5 zdtE^-^'^ 

AjJLE 




SxStHx^^ 


= -SxSt 
M 


A 

Uz^y 


-SyE^+--' 


1 (^^)^A AAA 
+ — i OxOzOxOtHz 
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6xEz^ ^ — SzEx 



SyEx'^^ — SxEy 



^ ^ -SxSxSyStHx'^^ 



4^£ 



SxSySzStHy^^ 



Removing the operator dx from all the terms in the above equations leads to the 
equivalent form of the ADI-FDTD scheme which is omitted here for shortness. 

From the above equivalent scheme and the Taylor expansion it can be seen that 
the truncation terms in (|4.2l) - (|4.13p are bounded if the exact solution of the Maxwell 
equations is suitably smooth, that is, 

(4.14) ISxP'lJ + |4/3«3 J + \SxPlJ < CehiiAtr + {Axf + {Ayf + {Azf}, 

(4.15) 14^,3, J + + \^^C,jJ < CeHiiAtf + {Axf + {Ayf + (Az)^}, 

where I ~l + ^ for I = k,n and C^h is a positive constant depending on and 
the upper bounds of the derivatives of e, h. 

For n > 0, w = x,y,z and for a = i, i + 1/2, /3 = j,j + 1/2, j — k,k + 1/2 let 



r cn A A P" 



where = (£■;', 5^, f^') and = {■H^,n'^,H'i). Then subtracting the 5^-ADI- 
FDTD scheme (i.e. the ADI-FDTD scheme applied by 6x) from the discrete form 
of the Maxwell's equations (see (|4.2p - (l4.13p ') leads to the following system of error 
equations: 
Error-(5^-Stage 1: 



Sx^: 



2e 



SxSyHz ^ — S^S!^ 



SxK 



At 
At 



dx£. 



y 



r _ At r c nj'r' 

Ox^y ^^^x^z'^x 



A ^"+1 _ A 0/"+^ — A 



SxSzC ' = SxK 



^^A A -OV" . ^*A 
-i^OxOztLy + —OxPx 



At 



-SxSxm 



At. 



2e 



y 



2e 

^SxSyS'!: + ^6xC-' 
2fi ^ ^ 2 



i-J,k 



2 2 '"•f^ 2 



At 



Ai e 



^.^r^ - T^'^.^c"^ - Sx-H'y ~ TT^xSzs: + 



2^ 

At 

2^ 



Sx'H-z ^ — — — 5x5y£: 



11+ i 



2^1 

SxK - ^SxSxS;; + ^SxC^' 



i,i+i;,k 



Error-(52;-Stage 2 can be obtained similarly. For example, the first error equation is 
given by 



2e ^ 2e 2 



We also need the following system of error equations for the ADI-FDTD scheme, 
which is obtained by repeating the argument above and letting Sx — I (the identity 
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operator): Error- Stage 1: 





2 




2 " 


2 




- 


2 




2 


2 






At „+i 


= e£^- 








MT^r^ - 




= t^n:- 




2 




mK^^ - 








yuAi „+l 
2 




/xT^r^ - 


2 "y'-'^ 






A* At „+i 
2 





Error-Stage 2 can be obtained similarly. For example, the first equation in this stage 

is given by 



In the above error equations, ^ , 'H™^ ^ and with w = x,y,z,m = n, h, 
can be regarded as those in Error-(5a;-Stages 1 and 2 with replaced with the identity 
operator /. For example. 



£•"-1-5 
Ox 



K 



U-t-^,J,fc' 



i,n* W"! 



n+1 

V 



e'y*, e"*, h"* and hy* can be similarly defined. Note that /SZ^^ and ^Z^^ with 
w = x,y, z are the truncation error terms for each equation of the equivalent form of 
the ADI-FDTD scheme (i.e. the equivalent form of the J^-ADI-FOTD scheme above 
with 5x — I)- For example, 



n+i. 



Ste. 







- Szhy^^ 


e 




U. 




- S e"+^" 

UyCz 





— OydxOtey 

Aij,e 



6zdx6th 



n+4 



If the solution of the Maxwell equations is suitably smooth then the following 
estimates hold 



(4.16) 



I + 1/3^" I + \^^A I < CeHiiAtf + {Axf + {Ayf + (A^)^ 



(4.17) |^:+',| + l^;;;:^,! + l^i^^^.l < CeuliAtf + [Axf + [Ayf + {Az)% 

where Ceh is a positive constant depending on £,/x and the upper bounds of the 
derivatives of e, h. 
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we now estimate (|4.1[) by making use of these error equations. Multiplying both 
sides of the first and second three equations in Error-J^j-Stage 1 by ^/e and respec- 
tively, and taking the square of the resulting equations, we will obtain six equations 
which are similar to the six equations (|3.10l) - (|3.15l) in the proof of Theorem 13.21 but 
each equation has one extra term /" (i = 1, 6) which related to the local truncation 
error term. For example, the first equation is given by 



where 



The other terms /f {i — 2, 6) can be defined similarly and obviously. 

Similarly, we can obtain six equations from Error-J^j-Stage 2. For example, the 
first equation in this stage can be obtained from (|4.16p : 



i+^,j,k 



where 

4 ze " 



The other terms f "^^ {i = 2, • • • ,6) corresponding to the other five equations are 
similar to /"^^ 

Using exactly the same argument as in the proof of Theorem 13.21 we arrive at 



(4.18) ll^r^+'lll + \\s,n''+'\\l + ^-j^{\\sj'^£''+'\\l + 

= ||<5.£"||| + \\6,nm + ^-P^iWsAml + IIMs'H"!!!) 

4:fie 

"II? + iiH-il? + ^ mn-\\j + \\6'^n\'i) 

+F"'"+^Ai; + G"'"+^AyAz, 
where Av = AxAyAz and for n > we have 

F-^' = E E E [/i - f^^\,. + E E E [/^' - f^''U..i. 

i=l j = l fe=l 1=1 j=0 k=l 

1-2 J-1 K-1 1-2 J-lK-l 

+ E E E [•^3" - z?^^] j,fe+i + E E E [•^4 - z"^^] 1 j+ 1 i 

1=1 j=l fc=0 i=l i=0 fc=0 
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/-1,7-l/f-l /-IJ-lA'-l 

'X! X! X! [■^5' ~ f5~^^'\i,],k+^ + X! X! X! Ife - f6~^^]i.j+} 

2=1 j = l fe=0 i=l j"=0 fc=l 



'n,n+l 



G 



Gn,n+1 
1 



2 +G-"-^)/Ax, 



i' = l,/-l j = l k=0 

E EEH^r+^;)- 

i' = l,/-l j=0 k=l 

GT^'= E EE 

i' = l J-l i=0 fc=0 



g; 



'n,n+l 



(At)! 
2 

(At)2 



,11+ J 



2 

(Ai)2 



n+1 7"^^ 



i' d.k 



,j,k 



i',j,k 



We now estimate each term in and (^"."+1. Noting that by the initial condition 

H° = £° = 0,we have 

n n /-1,7-1A'-1 

E^i = EEE E eAt[5A£'^^ + £i) + ^uMr-K) 

1=0 ;=o j=i j=i fe=i 
7-1 .7-1 Tf-l 



^ E E E [^^/^(^^-^^^r^)' + ^('^^^^'^^ 

i=l j=l fc=l 



+Ai5;]{e(<5,£:^) 



(=0 



4^£ 



i,j,k 



(4.19) +G{(At)4 + (Aa;)4 + (Ayf + (Az)*}. 



Similar estimates can be obtained for /" with i = 2, 6. Add these estimates together 
and use the definition of the discrete norms we get 

1 (A<)2 



1=0 



yji+l ||2 



{\\S,£-+'\\l + \\6,W-' ,H) 



(4.20) +C{{At)* + (Ax)* + {Ay)^ + (Az)^}. 
Arguing similarly as in deriving (|4.19p we can derive that 

(4.21) ±Gr'< E EE[^(^M(«"^)^ 

/=0 i' = l,7-l j=l k=Q ^ 



1=0 



4^e 



i' ■J,k+i 



C{A)\ 
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J-1 K-1 



(4.22) ^4''+^ < E E E [^(^M'5.Hr^)^+s(fr^)^) 

C(A)4; 



i=0 



i' = l,/-l j=0 /c=l 
n 

+AtE((4)' 



(At)2 



4/ie 



(4.23) EG3'+'< E EE[^(^^(^.^"^^)'+/«^^ 



n 

-At^(MHi)2 



(=0 



4^e 

4^£ 



i'j+5,fe+l 



+c•(A)^ 



where (A)"* = (At)* + (Aa;)"* + (Ay)"' + (Az)''. Combining (IT^ - ip:^ and using the 
definition of the discrete norms we obtain that 



1| 2'j 



ri+1 



115} 



(4.24) 



4/ie 



By summing up both sides of (|4.18p over the time levels n and using (|4.20p and 
(liTM)) it fohows that 

(4.25) + + M!(||4^J^£:n+i||2^ + P.^a^H^+illl) 



4/ie 



(At) 



2 r 



4/ie 



i=0 



4/ie 



'III, + 11452^-W'lll + ll^^'^'lll + \\5Wi) + C(A) 



The discrete Gronwall's Lemma 14.11 implies the estimate (I4.ip with w = a;. This 
completes the proof. □ 

If we apply the difference operators 5t and 5w5t with w — x^y, z to the equations 
in the ADI-FDTD scheme and argue similarly as in the proof of Theorem 14.21 then 
we are able to obtain the following result. 

Theorem 4.3. Suppose the exact solution e, h satisfies the same conditions as 
in Theorem \4-S\ Then for any n > and any fixed T > there is a constant C 
independent of At, Ax, Ay, Az such that for any w = x,y, z, 

ll^..5tf"+5||| + \\8^btU^^H\ + M!(||55.^^5^£"+i||2^ + ii4'^-'5tH"+^lll) 



|.o,n+i||2 ,(At)2 



l'52''^*H"+^||i(„) 



|(5f<5t£:"+=| 



Liw) 



(4.26) < C{{Atf + {Axf + {Ayf + (Az)*}, 

where L{x) — I, L{y) — J and L{z) = K . 
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Proof. We only prove (I4.26|) for the case with w — x. The other cases with w — y,z 
can be proved similarly. 

Repeating the argument in the proof of Theorem 14.21 with n being replaced by 
n — 1/2 and each function U being replaced by 5tU, we are able to obtain that for 
any n > 1, 

+ \\5AW'+--fH + ^(||<5.<5i"<5tf"+^||2^ + 1144' III) 
+||<5i£"+^||? + IIW+^II? + ^{||4W+^||? + \\5'l5tE^+'^\] 



+\\sts^ri + im^ri + ^{ii'^t^^^ii? + ii'^i'^t^^ii? 

+ (At)4 + (Ax)^ + (Ay)* + (Az)4}. 
Since = n° = 0, 6t£^^^ = {l/At)£'^ and St'H'^/^ = {1/At)n\ it is enough to prove 

(4.27) ||<5,f illl + WS^HYh + ^-^{WSAsYh + ll^^^'lll) 

< C{At f [{At)^ + {Ax)^ + (Ay)4 + (Az)^] . 
Letting n = in (|4.18p and using the fact that £^ = TiP = we find that 

(4.28) ii^s^lll + Ws^n^Wl +^-^{\\5A£'\\l + ll'^.'^^'^'lll) 



l^^ll?+IIH^II?+^(ll4'H^II?+p?^^ll? 



< F^'^Av + G'^'^AyAz, 

where Av — AxAyAz. We now estimate each term of and G"'^ slightly differently 
from the proof of Theorem 14.21 (cf. (j4.19p - (j4.24p ) since we do not need to take the 
sum with respect to n here. For example, taking n = in ()4.19p we have, on noting 
that £° = 0, that 



1=1 j=i k=i 



^EEE[^^MM.H,j)^ + f(4£^)^ 



-2 4u£ ■ ■ 2' U,j.k 

1=1 j=i fe=i ^ 

+C{At f [(At)4 + (Ax)4 + {Ay)* + (Az)*] . 

The other five terms /°, i = 2,..., 6, of F^'i and G-'\ i = 1,2,3, can be estimated 
similarly. Then by the definition of the discrete norms it is derived (cf. (I4.24p ) that 

F'''Av + G^''AyAz<l^-^{\\Sj'^£YH + \\Sj!^nTE 
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+\\sl£'\\j + \\5'^n'\\j} + -{\\6,£Te + Ws.nTH + W^'W! + \mj) 

+C{At f [{At)^ + (Ax)4 + {Ay)^ + (Az)*] . 

This together with (|4.28p implies the estimate (|4.27p . The proof is thus complete. □ 
Remark 4.4. It is a superconvergence result that the ADI-FDTD scheme is 

second-order convergent in space under the discrete semi- norm (see Theorems 14.21 

and l4.3|) . This is consistent with the superconvergence result of the semi-discrete Yee 

scheme established in [18] since both the ADI-FDTD scheme and the Yee scheme are 

based on the same spatial discretization technique. 

Setting Sw ^ I with w = x,y,z in Theorems 14.21 and 14.31 gives the following 

optimal error estimates in the discrete norm, which can be proved by arguing 

similarly as in the proof of Theorems 14.21 and 14.31 with Siu = /. 

Theorem 4.5. Suppose the exact solution e, h satisfies the same conditions as 

in Theorem \4-^ Then for any n > and any fixed T > there is a constant C 

independent o/ At, Ax, Ay, Az such that 

\\£-^'\\% + wh-^Th + ^-^(WS'^s'-^Yh + \\s'^n-+'\\l) 

(4.29) < C{(At)4 + {Axf + {Ay)^ + {Azf}, 

(4.30) < C{{At)^ + {Axf + {Ayf + {Azf}. 



5. Divergence preserving property of the ADI-FDTD scheme. As a con- 
sequence of the Maxwell equations (|2.2p - (|2.4p . if the electric field e and the magnetic 
field h (multiplied with e and /i respectively) start out divergence free, they will re- 
main so during wave propagation at any time and any place in the domain £7. In this 
section, we will prove that this property is preserved with the second-order accuracy 
in both space and time by the ADI-FDTD scheme. 

Theorem 5.1. Let E" and H" he the solution of the ADTFDTD scheme and let 
V'' • (eE") and V'' • (/iH") denote their discrete divergence: 

• (eE") = e{5,E^ + SyE^ + 5.i?,")|,^^. 
V'' • (a^H«) = ^l{^,H:^ + SyH^ + S.H^)\,. -^. 

Then 

(5.1) II V'' • (eE")||2 < C{||V'' • (ee°)||2 + (At)^ + (Ax)^ + {Ayf + (Az)*}, 

(5.2) II V'' • (mH")||2 < C{||V'' • (/ih°)||2 + {Atf + (Ax)^ + {Ayf + (Az)^}, 

/-I ,7-1 K-1 

where for a grid function U, \\U\\^ ~ elUi.j^k)"^ Av. 

1=1 j=l k=l 

Proof. We only prove (|5.1I) . The inequality (I5.2p can be proved similarly. 
From the property that V • (ee) = V • (ee°) it follows that 

II V'^ • (ee")||2 < ||V'^ • (£e")||2 + C{{Axf + {Ayf + (Az)^}. 
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Thus, and by Theorem 14.21 it is derived that 

IjV'' . (eE")f IIV'* • [e(e" - E" + e")]|p 

<2(||V''-(e5")||2 + ||V''-(£e")||2) 

< C{||V'' • (ee")||2 + {At)" + {Axf + {Ay)" + {Az)"} 

for some generic constant C > 0, which is (|5.ip . The proof is thus complete. □ 



6. Numerical experiments. In this section we carry out numerical experi- 
ments to verify the energy identities, optimal error estimates and divergence preserv- 
ing properties of the ADI-FDTD scheme. 

6.1. An exact solution and its energy conservation. Consider the Maxwell 
equations (|2.2p - (|2.4p with e = fi — 1. The exact solution is given by 



Sx = ^ cos(A/37ri) cos[7r(l — x)] sin[7r(l — y)] sin[7r(l — z)], 

= -VS cos{V3TTt) sin[7r(l — x)] cos[7r(l — y)] sin[7r(l — z)], 

Cz — — -VS cos{V3T:t) sin[7r(l — x)] sin[7r(l — y)] cos[7r(l — z)], 

hx = sin(-\/37rt) sin[7r(l — x)] cos[7r(l — y)] cos[7r(l — z)], 
hy — sin(V37rt) cos[7r(l — x)] sin[7r(l — y)] cos[7r(l — z)], 
hz — - sm{ViTTt) cos[7r(l — x)] cos[7r(l — y)] sin[7r(l — z)]. 

It is easy to verify that the above solution has the following energy conservation 
and divergence-free properties: 

llelp = ^cos2(V37rt), \\h\f ^ ^ sin^ {VSnt) , V • e = V • h = 0, 
64 64 

p2 1,5^ 2 ah 2 21 2 p2 II 2 , II 2 68 4 

where for a vector function u {u^, Uy,Uz) the L^-norm || • || is defined by 

||u||^ = ^ P{n)(\ux\^ + \uy\'' + \uA^)dx, P{e) - e, P{e) = ^. 

Eehw and Eehwt can be defined similarly. Note that Eehw = E^hx, Eehwt = Eehxt, 
w = y,z. 

6.2. Energy identities of the ADI-FDTD scheme. Denote by | • | i and 
II • II 2 the two new norms defined in Theorems 13.21 and 13.41 respectivelv. Then 

||(E",H")||? = ||(5xE"||| + ||4H"||lf + ||E"||? + ||H"||? 
, (At)2 



■ {WSxS'^E'^rH + ll'5.<52'H"||| + ||52'H"||? + ||<5i^E"||2 
(E",H")||2 ^ ||E"||| + IIH-lll, M!(||j^^H"||| + ||^^E"|||, 
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Similarly, we can define ||((5tE", 5tH")||2 and ||(5fE",5tH«)||2. To verify the four 
energy identities established in Subsection 13.31 we introduce the following notations: 

||(E",H")||,-||(E°,H°)||, IKE",H")||i IKE'\H")|b 

||(EO,HO)||j ' """"i E,u. ' ' ' 

||(^,E"+i<^,H"+^)||,-||(,5,EiH^)|l, ^ ||(^,E^^,H'^)||, ^ 

||(5,Ei,5,Hi)|[, ' ' ■ 

Here, EH^^'^, EHl'j*, j — 1,2, represent the relative error of the two new energy norms 
of the solution to the ADI-FDTD scheme at time level i" and at the initial time or 
t^, and EHj^, EHJ^p j = 1, 2, stand for the ratio of the two new energy norms of the 
solution of the ADI-FDTD scheme to the corresponding energy norms of the exact 
solution. 

The computational results are presented in Table 16.11 for different time levels n 
(n = T/At), where the solution to the ADI-FDTD scheme is computed with the 
spatial step sizes Aa; = Ay = Az = 0.01 and the time step size At = 0.01 (note that 
the Courant number is -\/3) • From Table 16.11 it is seen that the four discrete energy 
identities shown in Theorems l3.2[ 13.31 13. 4l and l3.5l are satisfied and the discrete energy 
of the approximate solution is almost equal to the corresponding ones of the exact 
solution (their ratio is almost equal to 1). 

Table 6.1 

Relative error of the discrete energy of the ADI-FDTD solution at time level n and at the initial 
time, and the ratio of the discrete energy of the ADI-FDTD solution to the true energy of the exact 
solution under the norms \\ ■ ||i and \\ ■ ||2 at different time levels with Ax = Ay = Az = At = 0.01 



Energy \ Time 


n=100 n=400 n=800 n=1600 n=2000 


EHf 
EH^ 


3.157e-13 4.780e-14 2.783e-13 2.537e-13 1.646e-13 
1.00008 1.00008 1.00008 1.00008 1.00008 


EHlf 


1.627e-12 1.518C-12 1.428e-12 1.478e-12 1.559e-12 
0.99979 0.99979 0.99979 0.99979 0.99979 


EH'^^ 
EH^ 


1.397e-13 2.808e-13 3.432e-13 2.254e-13 2.167e-13 
1.0001 1.0001 1.0001 1.0001 1.0001 


EHlf 
EHl\ 


1.457e-12 1.585e-12 1.744e-12 1.438e-12 1.634e-12 
0.99984 0.99984 0.99984 0.99984 0.99984 



EHf = 



Em 



To verify the energy identities without the perturbation terms, define 

||(E",H")||?, := ||<5,E"||| + ||<5,H"||2, + ||E"||2 + ||H"||2, 
|!(E",H")||i :-||E"||| + ||H"||l,. 

Similarly, let EH^f, EH^jl, j = 1, 2, represent the relative error of the two new energy 
norms without the perturbation terms of the solution to the ADI-FDTD scheme at 
time level t" and at the initial time i° or and let EHj\, EH^j^, j = 1, 2, stand for 
the ratio of the two new energy norms without the perturbation terms of the solution 
of the ADI-FDTD scheme to the exact ones. 

Table \6l2\ shows similar results to those in Table IHTTl with the four discrete energy 
identities without the perturbation terms. These new discrete energy norms can be 
seen as the discrete forms of the physical energies. From these values we can see that 
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the energy identities without the perturbation terms are still satisfied and that the 
discrete energies of the ADI-FDTD solution are almost equal to the corresponding 
ones of the exact solution. 

Table 6.2 

Relative error of the discrete energy of the ADI-FDTD solution at time level n and at the initial 
time, and the ratio of the discrete energy of the ADI-FDTD solution to the physical energy of the 
exact solution under the norms \\ ■ and || ■ ||2* at different time levels uiith Ax = Ay = Az = 
At = 0.01 



Energy\Time 


n=100 n=400 n=800 n=1600 n=2000 


EH[\ 


1.592e-14 3.140e-13 1.730e-14 1.413e-13 1.070e-14 
0.99996 0.99996 0.99996 0.99996 0.99996 


EHl\% 


1.942e-13 3.359e-13 2.508e-13 2.899e-13 4.480e-13 
0.99967 0.99967 0.99967 0.99967 0.99967 




1.463e-13 1.421e-13 1.531e-13 2.105e-13 6.590e-15 
1.0000 1.0000 1.0000 1.0000 1.0000 


EHlf, 
EH^2. 


9.150e-14 9.050e-14 2.309e-13 2.269e-13 3.705e-15 
0.99971 0.99984 0.99971 0.99971 0.99971 



6.3. Error and convergence rate of the ADI-FDTD scheme. To verify 
the convergence rate of the ADI-FDTD scheme under different norms, let 

E-H'i = ||(£",H")||i/£;e/.., EU^^ = \\{5tE^+^-,6tW'+H\i/Eeh.u 
E-H'i = ||(£",H")||2/£;eh, EWi^ = mE'^+'-^StU^+HU/E^hu 
EK^^{\\E-\\l + \ml)/E,u. 

Table 16.31 gives the relative error in different norms of the ADI-FDTD solution 
at different time levels with the time step size Ai = 0.01, the spatial step sizes 
Aa; = Ay = Az = 0.01 and the Courant number v^- The results demonstrate that 
even for a large Courant number and at a long time, e.g. T = 20, the ADI-FDTD 
solution is still stable and accurate in different norms. 

Table 6.3 

Relative error in the norms \\ ■ ||i, and \\ ■ ||2 of the ADI-FDTD solution at different time with 
Ax = Ay = Az = 0.01 and At = 0.01 



Error\Time 


n=100 n=400 n=800 n==1600 n=2000 


EH1 


9.091e-4 3.580e-3 7.161e-3 1.432e-2 1.790e-2 
9.162e-4 3.585e-3 7.159e-3 1.431e-2 1.789e-2 


En^ 

^2 


3.214e-4 1.266e-3 2.532e-3 5.063e-3 6.329e-3 
9.159e-4 3.585e-3 7.195e-3 1.431e-2 1.789e-2 



The convergence rates in time and space of the ADI-FDTD scheme are given in 
Tables 16.4116.91 under different norms. The results show that the convergence rate in 
the norms || • ||i, || • ||2 and || • ||o is second order in space and time, which confirm the 
theoretical results. Note that the second-order convergence in time and space under 
the norm || • ||i is a superconvergence result. 
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Table 6.4 

Convergence rate in time under the norm || ■ ||i at T = I with Ax = Ay = Az = 0.01 



At 


Sn^i Rate 


Sn^i Rate 


0.05 
0.04 
0.025 
0.02 


1.749e-2 
1.127C-2 1.971 
2.978e-3 1.938 
2.978e-3 1.881 


1.706e-2 

1.108e-2 1.934 
2.975e-3 1.913 
2.975e-3 1.863 



Table 6.5 

Convergence rate in space under the norm || ■ ||i at T = 1 with At = 0.001 



Ax = Ay = Az 


£-Hf Rate 


Rate 


0.025 
0.02 
0.01 


1.405e-3 

9.016e-4 1.987 
2.305e-4 1.968 


1.430e-3 

9.221e-4 1.966 
3.141e-4 1.554 



Table 6.6 

Convergence rate in time under the norm \\ ■ \\2 at T = I with Ax = Ay = Az = 0.01 



A< 


Sn^ Rate 


£n]^2 Rate 


0.050 
0.040 
0.025 
0.020 


6.185e-3 
3.984e-3 1.971 
2.978e-3 1.938 
1.053e-3 1.881 


1.706e-2 

1.108e-2 1.934 
2.975e-3 1.913 
2.975e-3 1.863 



Table 6.7 

Convergence rate in space under the norm \\ -112 at T = I with At = 0.001 



Ax = Ay = Az 


Sn'^ Rate 


Sn'^^ Rate 


0.025 
0.020 
0.010 


4.968e-4 
3.188e-4 1.988 
8.149e-5 1.968 


1.428e-3 

9.165e-4 1.987 
2.342e-4 1.968 



6.4. Divergence of the numerical electric field. To verify the discrete 
divergence-free property of the ADI-FDTD scheme, we introduce two norms for di- 
vergence: 

Divioo 
Divi2 

Table I^TTUl gives the values of the above norms at different times T — 1,4,8, 16 and 
20, where MAt = T. The results show that the discrete divergence is almost zero, 
which means that the ADI-FDTD scheme is of the approximate divergence preserving 
property. 
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Table 6.8 

Convergence rate in time under the norm \\ - Wo at T = 1 with Ax = Ay = Az = 0.01 



At 


Rate 


||-H"|l£; Rate 


0.050 
0.040 
0.025 
0.020 


1.914e-2 

1.232e-2 1.975 

4.955e-3 1.938 
3.260e-3 1.876 


1.595e-2 

1.031e-2 1.956 

4.157e-3 1.933 
2.730e-3 1.883 



Table 6.9 

Convergence rate in space under the norm \\ ■ ||o at T = 1 with At = 0.001 



Ax = Ay ^ Az 


WS'^'We Rate 


WH'^We Rate 


0.025 
0.020 
0.010 


4.968e-4 
3.188e-4 1.988 
8.149e-5 1.968 


1.7064e-2 
1.1082e-2 1.934 
2.975e-3 1.913 



Table 6.10 

Divergence error of the numerical electric field E in L°° and I? norms at different times 



Error/Time 


T=l 


T=4 


T=8 


T=16 


T=20 


Divioo 


7.060e-13 


1.421e-12 


1.943e-12 


2.956e-12 


2.912e-12 


Divi,2 


7.421e-14 


1.540e-13 


2.176e-13 


3.085e-13 


3.451e-13 
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